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Abstract 

Long-wavelength gravitational waves can induce significant temperature 
anisotropy in the cosmic microwave background. Distinguishing this from 
anisotropy induced by energy density fluctuations is critical for testing infla- 
tionary cosmology and theories of large-scale structure formation. We describe 
full radiative transport calculations of the two contributions and show that 
they differ dramatically at angular scales below a few degrees. We show how 
anisotropy experiments probing large- and small-angular scales can combine 

to distinguish the imprint due to gravitational waves. 
PACS NOs: 98.80.Cq, 98.80.Es, 98.70.Vc 
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The cosmic microwave background (CMB) temperature anisotropy may be induced by 
energy density fluctuations and by gravitational waves jl| , corresponding to scalar and ten- 
sor metric perturbations, respectively. Although anisotropy measurements probing angular 
scales above a few degrees (e.g., COBE 0) cannot discriminate scalar from tensor ||, we 
show in this Letter that the two contributions can be separated when data from smaller 
angle experiments are used as well. 

Resolving the two contributions relies upon detailed theoretical predictions for the form 

(S) (T) 

of the multipole components, a) m and a tm , of the relative temperature pattern on the 
sky, AT/T(9,4>). For inflationary models, each multipole for the two modes is predicted 
to be statistically independent and Gaussian-distributed, fully specified by angular power 
spectra, Cj, S ^ = (\af2\ 2 ) and C% = {\a^\ 2 }. Although Cf^ has been calculated before 
Q, was previously known only for low multipoles, £ ^ 30, relevant for angles above 
a few degrees |l]J^,§]. For these ts, the dominant spatial wavelengths contributing to C\ 
and C\ were outside the horizon at photon decoupling, and both scalar and tensor modes 
induce similar red shifts and blue shifts in the CMB [0,^,0. For example, COBE's DMR is 
unable to distinguish the two contributions due to cosmic variance (from the theory signal) 
and experimental noise, although their sum can be determined from one year of DMR data 
to within 30% (see Fig. 2a below), improving to about 15% with four full years of data. 

(T) 

In this Letter, we compute C\ to much higher multipoles, and show that the predicted 

(T) (S) 

C\ becomes highly suppressed relative to C\ at large i. The dominant wavelengths for 
I > 30 were inside the horizon at decoupling. Inside the horizon, scalar-mode anisotropies 
are enhanced by the gravitational instability of density perturbations and by Thomson 
scattering from moving electrons, whereas gravitational waves disperse as freely propagating, 
massless excitations and red shift away. Taking advantage of this difference, we find that 
combining experiments at small and large angular scales can determine the scalar and tensor 
components. Current CMB anisotropy data at small scales is not yet good enough to do so, 
especially since some of the signals observed may be Galactic rather than cosmic in origin, 
but the statistical errors of these and other near-future experiments are small enough to 
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allow separation at about the two sigma level, as we show here. 

Separating tensor from scalar is essential for theories of cosmic structure formation since 

(S) 

it is from C\ that we can infer the amplitude of the primordial density fluctuations. It 
also provides a critical test for inflationary cosmology [||,|]||. All inflation models produce 
a post-inflation spectrum of scalar and tensor metric fluctuations with some tilt away from 
the scale- invariant (spectral index n s = 1) form. The degree of tilt depends upon the details 
of the inflation model but may be quite significant HHH (n s <: 0.5). For a wide range 
of inflation models, but not all, the gravitational wave content of the CMB anisotropy is 
related to the tilt ||.|U : 

CP 7(1 -n.) , (1) 

where C% and C% are the £ = 2 (quadrupole) components of the power spectrum. Con- 
firmation of Eq.(l) would provide detailed information about the kind of inflation that took 
place and, consequently, about the evolution of the very early Universe. 

For a given sum of and fixed by large angular scale measurements, a larger 
tensor component reduces the small-angle anisotropy. Several other factors can have a 
qualitatively similar effect {e.g., increased tilt, 1 — n s ; decreased baryon density, Qb] and 
a nonstandard recombination history), but there are quantitative differences which we now 
describe. 

To compute Cf\ we evolve the distribution function, f(x, q, t), for photons at posi- 
tion x at time t with momentum q, using first-order perturbation theory of the general 
relativistic Boltzmann equation for radiative transfer |3| , with a Thomson scattering source 
term. Photon polarization is included by making f a 4-dimensional vector with components 
related to the Stokes parameters (f s with s = t,p,u,v correspond to the usual I,Q,U,V 
Stokes notation) and applying Chandrasekhar's development of the scattering source term 
for Rayleigh (and thus Thomson) scattering in a plane parallel atmosphere [7[ . In the scalar 
case, only f t and the 'polarization' f p are needed, so 2 transfer equations are required ||. In 
the tensor case f u also does not vanish, but it is related to f p , so again only two perturbed 
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transfer equations turn out to be required. To describe these equations, we introduce the 
relative perturbed distribution functions || = A5f s /(T Q df/dT ), where T is the CMB 
temperature and / is the unperturbed Planck distribution. 

To evolve the coupled equations, both A^ T ) and the metric are expanded in plane waves. 
In the frame in which the wavevector k is along the z-axis, the gravitational wave degrees 
of freedom in the metric are the transverse traceless modes, h + = h\\ = — /122 and h x = 
hi2 = ft<2i, which obey an Einstein equation, the wave equation for free massless particles: 
h € + 2^h e + k 2 h € = 0, e = +, x. Here the dot denotes derivative wrt conformal time, 
t = J dt/a(t) where a is the expansion factor, solved for by evolving the Friedmann equation. 

The radiative transfer equations for the two gravity wave polarizations separate, having 
an overall factor of cos(20) for e = + and of sin(20) for e = x, where (6,<p) are the polar 
angles, which we remove by introducing new variables, following Polnarev || 

A (T) = A^ ) (l- / i 2 )cos20 + A^ ) (l-/i 2 )sin20 

ACT) = Ag>(l + li 2 ) cos 20 + Ag } (l + /1 2 ) sin20 (2) 
Af) = -A| t + ) 2/isin20 + A^ ) 2 / ucos20 

The combination A^) + A^p is unexcited by gravity waves, as is Af), so the 4 Stokes 
radiative transfer equations reduce to two: 

A^P = —ikfiAfr — h e — a<7T n eA-te + acr T n e \l/ e , 
Ape = —ikfiA^p - aa T n e A^ - aa T n e ^ e , 

1 A ( T ) 1 A (T) 1 3 A ( T ) ^ 



IQ^tefi 7^6,2 + 70^,4 

3a( T ) 6a( T ) 3 \(T) 
5^V,0 7^6,2 70^pe,4 

Here <jt is the Thomson scattering cross- sect ion, and n e , the free electron density, is evolved 
using a careful treatment of the recombination atomic physics. 

As in the scalar case [|]], we solve these equations by expanding in Legendre polynomials, 
e.g., A te = £^(2£ + l)A te /Pe(Li), converting (3) to a hierarchy of coupled equations. Our 
solutions, A t e/(k, r) and A pt ^{k, r), can be combined into the power spectrum by summing 
over k and polarizations: 



C { P = f £ k (l - 1)1(1 + !)(* + 2)([ pt^li) 



(4) 



'(2^-l)(2£+3) 1 (2^+l)(2£+3) 

where AiG,£ = (^t+,e — iA tx ,e)/ A similar expression applies for the polarization power 
spectrum. || A useful check is to assume recombination is sudden at r = r r . The free- 
streaming solution from r r to the present r is then A t€ ^ = (—i) e drji(k(TQ — r))h e (r). 

( r p\ 

Substitution into eq.(4) gives the Abbott and Wise [0] formula for C\ , with which our 
numerical results agree for low £. 

In this paper, we discuss results for standard cold dark matter (CDM) models with 
normal recombination, although it is straightforward to adapt the numerical codes to other 
cosmo logical models (e.g., mixed hot and cold dark matter). We let n s , the ratio of tensor- 



to-scalar quadrupole anisotropy, the baryon density Qb and the Hubble parameter |10[ h 
freely vary. 

Fig. 1 shows a CDM model with Qb — 0.05, h = 0.5 and n s = 0.85, which accounts for the 
slight downward tilt in C\ at small I. (We plot £(£ + l)Ce, since it is flat for scale-invariant 
(n s = 1) Cf S ^ at small £.) The sharp increase in for £ > 50 followed by increasingly 
damped oscillations are due to adiabatic compression of photons and Doppler shifts during 
decoupling [|]]. The tensor mode behaves quite differently. The first moments drop sharply; 
then the curve settles to a tilted spectrum similar to the scalar case for 5 > £ > 50. At 

(T) ( 

£ <: 50, the tensor drops sharply just as the scalar rises. We have set Of 7 Of « 1, the 
inflation prediction of Eq. 1 for a tilt of n s = 0.85. We also illustrate how a scalar-only 
spectrum with low Qb (e.g., the dashed curve) partially mimics the scalar-plus-tensor shape 
for Qb = 0.05, (assuming a larger C\ )• Clearly, precise measurements are required to 
separately determine C^' jC^' and Qb- 
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FIGURES 

FIG. 1. (a) Angular Power Spectra for the tilted standard CDM model shown for tensor, 
scalar and the sum. The light dashed line is an Ob = 0.01 model, (b) shows the filters for the 
experiments used in this paper as examples (heavy lines). The light lines are other representative 
experiments. 

The filter functions, which indicate the experimental sensitivity as a function of t, are 
shown for various experiments in Fig. lb. The theoretical prediction for the rms fluctuations 
at each of the points in an experiment is found by multiplying Ce/4ir by the filter and 
summing over t, m. From Fig.l, we see that large-angle experiments (e.g., DMR ||, MIT [1 1 



and Tenerife |]12|) are equally sensitive to tensor and scalar modes, smaller angle experiments 



e.g. SP89 [13| and OVRO [16]) are sensitive mostly to scalar, while the intermediate SP91 



[14], 15] can measure some tensor, although predominantly scalar. 



A quantitative experimental fit to cosmological parameters is obtained by construct- 
ing likelihood functions C e (C2 > C2 , n s,^B,^) for each experiment, e, assuming Gaussian 



statistics. [18] Assuming the experiments are statistically independent (because they cover 
unrelated regions of the sky or very different angular wavebands), we combine the £ e 's to get 
the full likelihood, £ = H e C e , as shown in Fig. 2. For all but DMR, £ e is calculated using 
Bayesian techniques (|HJ] which take into account the removal of any linear combinations 
of the data such as gradients or averages by marginalizing over the coefficients, assuming 
uniform prior probability distribution in these coefficients. For SP91, the method was ex- 
tended to treat simultaneously the 4 frequency channels. For DMR, we used the Smoot et 
al. '90 A+B X 53 A+B' (quadrupole-subtracted) correlation function with a Gaussian 
approximation for the likelihood |19| A more complete analysis will only be possible once 
the DMR data are released. In Fig. 2, we have taken VLb = 0.05(2h)~ 2 , consistent with 
nucleosynthesis limits pCfl , and h = 0.5. 

Fig. 2(a) displays the DMR likelihood contours in the C 2 C\ plane for fixed n s 
(0.85), demonstrating that DMR can measure C2 + C2 , but cannot discriminate scalar 
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and tensor. A preferred tensor-scalar ratio does arise in Fig. 2(b) as soon as we incorporate 
small angle data. Figure 2(b) combines the DMR data, the 4- frequency-channel data from a 
9 point strip |14j] and a 13 point strip [|15|] in the SP91 experiment, as well as SP89 [13] and 



OVRO |T6[ data (which give weak upper limits but no detections). The SP91 scans appear 
to have detections, but the signal may be contaminated by unknown sources. 

Fig. 2(b) is tantalizing but inconclusive evidence for a gravitational wave contribution. 
Future refinement can be anticipated by using simulated data sets, constructed by taking 
single realizations of theoretical signals and adding experimental noise associated with sta- 
tistical errors (but no systematic errors). In Figure 2(c) and 2(d), the input signal is for 
a standard CDM model with n s = 0.85 and equal Cf° and C^ S) ([7.5 x 1(T 6 ] 2 ). We then 
simulate a suite of plausible near-future experiments: DMR with 4-year error bars; six 13 
point strips from an SP91 configuration (\8-27fiK error bars for each of the 4 frequency 
channels); six 9 point strips from an SP89 |13[ configuration; and an OVR022 configuration 
(7' beam, 22' double-difference throw, with 25 fiK error bars). Reduced (~ 15fiK) error bars 
were taken for SP89 to represent ongoing or planned experiments with beams ~ 0.5° which, 
with multifrequency observations, can achieve these sensitivities [|T7| . 



Fig. 2(d) shows two projections onto the C^-n s plane. The heavy contours are the 
likelihood if jC\ is restricted to the trajectories predicted by inflation, Eq. (1). The 
maximum lies within 10% of the input signal. In contrast, the light curves show the contours 
when the (unrestricted) maximum likelihood value for the given Cf } and n s is taken. The 
extended 1 sigma band along C% /C^p ~ (2 — 5.2(1 — n s )) 2 indicates an inability to dis- 
tinguish large tilt from large tensor component for these sensitivities. The band runs across 
the inflation prediction (heavy lines), intersecting in a narrow range about n s « 0.83, very 
close to the input value. 

We conclude that current and near-future anisotropy experiments are unable by them- 
selves to definitively test inflation (Eq. 1) or determine the gravitational wave contribution 
to the CMB. For the short-term, conclusions can only be drawn by adding extra assump- 
tions and/or other data. For example, we have already shown in Figs. 2(c,d) that, if Eq. 1 is 



assumed, n s and the gravitational wave imprint can be determined to within 2 sigma. Alter- 
natively, other cosmological constraints can be invoked. For example, a variety of arguments 
imply that the rms amplitude of the density fluctuations on scales of 8h~ 1 Mpc (erg), which 
is used to measure the amount of nonlinear dynamics in large-scale structure calculations, 
cannot lie outside of the range 0.45 and 1 for CDM models j2"T,|5|; this translates to the 
shaded region in Fig. 2(d). If one assumes a model of galaxy clustering with linear biasing, 
a tilt n s < 0.65 is required for standard CDM models to reproduce the galaxy correlation 
function || (which would exclude our input model). This severe restriction may be relaxed 
with less simplistic CDM models of galaxy formation, and in cosmologies with more power 



in the density fluctuations than CDM has on large scales [21 



The long-term future is brighter. Extensive mapping of the microwave sky on small 
and intermediate angular scales can lead to highly accurate determinations of the spectrum; 
e.g., using the filter functions of Fig. 1, we find that the limiting cosmic variance uncertainty 
in AT/T is only a few per cent for SP89 and SP91 configurations. Even at large angles 
where the cosmic variance is higher, we find that 5% accuracy should be achievable. Hence, 
given optimal experimental designs for measuring large and small angle anisotropy, there 
should be sufficient resolution for a fully independent test for inflation, theories of large-scale 
structure, and the imprint of gravitational waves. 
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FIG. 2. Likelihood contour maps for scalar 
{[C { 2 S) ] 1 / 2 /1CT 5 ) v.s. tensor ([C,j, T) ] 1/2 /1CT 5 ) amplitudes are shown in (a),(b) and (c) for the stan- 
dard CDM model with fixed n s = 0.85 tilt. The light curves are 1, 2 and 3 sigma lines, the heavy 
curve or 'x' gives the maximum likelihood, (a) DMR only, (b) DMR plus the 9 and 13 point 
SP91 data, along with SP89 plus OVRO. (c) shows the maps with simulated large and small angle 
data consisting of DMR (with 4 year error bars), six 13 point SP91 strips, six 9 point SP89 strips 
and one OVR022 strip. The mean CDM signal input into the simulated data is denoted by the 
square, (d) Shows 1,2 and 3 sigma likelihood contours for the simulated data in [C^] 1 / 2 — n s 

(T) (S) 

space, constrained to the C 2 /C 2 trajectory defined by Eq.(l) (solid) and the unconstrained 
maximum likelihood trajectory (dashed). Shading indicates the range for which CDM models are 
not dynamically viable. 
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